Complete characterization of the macroscopic 
deformations of periodic unimode metamaterials of 

rigid bars and pivots. 

Graeme Walter Milton 

Department of Mathematics, University of Utah, Salt Lake City UT 84112, USA 

Abstract 

A complete characterization is given of the possible macroscopic deformations 
of periodic nonlinear affine unimode metamaterials constructed from rigid bars and 
pivots. The materials are affine in the sense that their macroscopic deformations 
can only be affine deformations: on a local level the deformation may vary from cell 
to cell. Unimode means that macroscopically the material can only deform along 
a one dimensional trajectory in the six dimensional space of invariants describing 
the deformation (excluding translations and rotations). We show by explicit con- 
struction that any continuous trajectory is realizable to an arbitrarily high degree 
of approximation provided at all points along the trajectory the geometry does not 
collapse to a lower dimensional one. In particular, we present two and three dimen- 
sional dilational materials having an arbitrarily large flexibility window. These are 
perfect auxetic materials for which a dilation is the only easy mode of deformation. 
They are free to dilate to arbitrarily large strain with zero bulk modulus. 

Keywords: A. microstructures, B. inhomogeneous material; B rods and cables; B foam 
material; unimode material 

1 Introduction 



Milton and Cherkaev (1995) (see also chapter 30 of Milton (2002)) introduced the con- 
cept of unimode, bimode, trimode, quadramode and pentamode materials in the context 
of linear elasticity. In these extremal three-dimensional materials the six eigenvalues of 
the elasticity matrix split into two groups: very large ones and very small ones: the 
number of very small ones (counting degeneracy) determines the category of extremal 
material: unimode if there is one small eigenvalue, bimode if there are two, trimode if 
there are three, quadramode if there are four, pentamode if there are five. For periodic 
pin-jointed trusses the category, and the infinitesimal modes of deformation, can be de- 



termined by an application of Bloch's theorem as noted by Hutchinson and Fleck (2006) 
An example of a unimode material is an elastically isotropic material with a Pois- 
son's ratio arbitrarily close to —1 so that the only easy mode of deformation is an in- 
finitesimal dilation. Isotropic materials with negative Poisson's ratio, the existence of 



1 



which was questioned for a long time, were first manufactured by Lakes (1987) Elas- 
tically isotropic materials with a Poisson's ratio arbitrarily close to —1 were found by 
lmgren (1985) and Milton (1992) The constructions of Almgren used sliding collars, 
whereas the constructions of Milton did not. There has since been extensive work on 
materials with negative Poisson's ratio which are known as Auxetic materials: see, for 



example Greaves et al. (2011) and Mitschke et al. (2011) and references therein. An ex- 



ample of an isotropic pentamode material is a gel: it is easy to deform under any five 
independent infinitesimal shears, yet strongly resists hydrostatic compression. In gen- 
eral, pentamode materials can resist a desired combination of compression and shear, 
and as such they are useful for transformation acoustics and in particular cloaking for 



acoustics, as discovered by Norris (2008) Designs for these more general pentamode ma- 
terials suggested by Milton and Cherkaev (1995) were recently experimentally realized 
by Kadic et al. (2012) dramatically illustrating the present day ability to tailor make 



designer microstructures (see also Biickmann et al. (2012)). 



The concepts of unimode, bimode, trimode, quadramode and pentamode materials 
can be extended to non-linear elasticity and find a natural place when one is interested 
in trying to find what metamaterials can be produced with two and three dimensional 
periodic arrays of rigid bars and pivots. Since they are periodic they have an underlying 
Bravais lattice. A two dimensional Bravais lattice consists of points 



x = iu + jv, 



as % and j range over all integers. Let F be the deformation matrix with the primitive 
vectors u and v as columns. Then under a rotation R, since u and v transform to Ru 
and Rv it follows that F transforms to RF leaving the symmetric Cauchy-Green matrix 
A = F T F invariant. There are many other Bravais lattices. In particular, there are 
lattices with lattice vectors 



u = ku + £v, v' = mu + nv, 



;i.2) 



for any choice of the integers k, £, m and n such that u' and v' are independent. The 
corresponding deformation matrix F', which has u' and v' as columns, is given by 



F' = FM, where M 



k m 
i n 



;i.3) 



and M is non-singular. Other Bravais lattices will correspond to matrices M with ele- 
ments that are not necessarily integers. As the material continuously deforms, with the 
Bravais lattice having lattice vectors u and v undergoing an affine transformation, F 



follows some trajectory F(t ) beginning at to 
with the deformation is a path 



t and ending at t = t$. Associated 



C(to) = [F(t c 



\Tl-l 



F(to) r F(t )[F(t )]- 



;i.4) 



in the three dimensional space of symmetric matrices beginning at C(£g ) — I- If we chose 
a different Bravais, such as that in (11.21) . this path, for the same deformation, would be 
the same since 



[M T F(t ) r ]- 1 M r F(t ) T F(t )M[F(t )M]^ = [F(t " ) T ]^F(t ) T F(t )[F(t )]~\ (1.5) 
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On the other hand there can be continuous deformations in which the larger lattice un- 
dergoes an affine transformation, but the smaller lattice is distorted. So it is possible 
that there exist additional paths of deformation C(io) for the larger lattice which are not 



accessible to the smaller lattice, see figure [T] and the examples in Kapko et al. (2009) 



for instance. If there exist deformations which are non-affine at the macroscopic scale, 
then the structure is said to be non-affme (see figure [U^c)). In this paper we restrict our 
attention to affine materials, i.e. materials for which the only macroscopic deformations 
are affine ones. More precisely a material is affine if and only if any deformation x'(x) 
(defined only for points x on the rigid bars which get moved to x'(x)) necessarily ap- 
proaches Bx + b as |x| — > oo, for some non-singular matrix B and vector b (dependent 
on the deformation) to within terms of o(|x|). A material which isn't affine is non-affine. 




Figure 1: A parallelogram array of rigid bars can deform from (a) to (b). The red dots 
denote the Bravais lattice which undergoes an affine transformation. As a result the array 
is not unimode, but instead trimode. It is a non-affine material as it has macroscopic 
deformations which are non-affine as in (c). 

A material is classed as null-mode if the only possible continuous deformation path 
for any Bravais lattice is the point C(to) = I. Some internal motions could be pos- 
sible but macroscopically the material is rigid to periodic deformations (in the sense 
that for any Bravais lattice the only continuous deformations that keep the periodic- 
ity of that lattice only rotate and translate the lattice). A material is classed as uni- 
mode if it is not null- mode and the possible deformations C(to) for any Bravais lattice 
all lie on the same one dimensional curve. A material is classed as bimode if it is 
not null-mode or unimode and the possible deformations C(t ) for any Bravais lattice 
all lie on the same two dimensional surface. Four examples of such two-dimensional 
unimode materials that may be constructed from bars and pivots are given in figure 
[21 with underlying Bravais lattices marked by the red dots. The kagome lattice also 
appears to be a unimode material for which the only macroscopic mode of deforma- 
tion is a dilation [see, for example, IGrima and Evans 2006( IHutchinson and Fleck 2 006; 
Kapko, Treacy, Thorpe, and Guest 2009 and Sun, Souslov, Mao, and Lubensky 2012] . 



Guest and Hutchinson (2003) investigated the question as to whether it is possible 
to build a periodic pin-jointed structure having no easy modes of deformation and such 
that replacing any bar by an actuator and changing its length leads only to a change of 
the geometry of the structure, and not to self stress. They pointed out this is of interest 
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to designing adaptive structures. They found that no such periodic rigid pin-jointed 
structures can exist. However the periodic unimode materials considered here (unlike 
the kagome lattice) have the property that they become macroscopically rigid if a single 
(appropriately placed) bar is added to the structure. Changing the length of this bar 
(within limits) leads only to a change of the geometry of the structure, and not to self 
stress. Replacing that bar by an actuator leads to interesting adaptive structures. 
A three dimensional Bravais lattice consists of points 

x = ill + j v + kw, (1.6) 

as i, j and k range over all integers, and we take F as the matrix with the primitive 
vectors u, v and w as columns. Then as the Bravais lattice undergoes a continuous 
affine transformation, the matrix F varies along a path F(io) beginning at to — an d 
ending at to —t^, and the symmetric matrix C(io) defined by (11.41) follows a trajectory 
in six dimensional space. In three dimensional unimode, bimode, trimode, quadramode 
and pentamode materials the possible trajectories C(t ) are confined to, respectively, a 
one, two, three, four, or five dimensional manifold, for all Bravais lattices and for all 
possible affine deformations of those lattices. 

Here we only consider unimode materials built from rigid bars and pivots and address 
the problem of characterizing what trajectories C(io) are realizable and in the process 
we find multiscale microstructures which realize given trajectories C(io) to an arbitrarily 
high degree of approximation. We do not consider trajectories where the microstructure 
degenerates into a lower dimensional geometry, and thus we assume C(to) remains strictly 
positive definite, and so the trajectory C(£o) is confined to lie within the "cone" of positive 
definite matrices. Here we essentially show, by construction, that every such trajectory is 
realizable to an arbitrarily high degree of approximation: there are no hidden constraints. 

Most of the structures we consider here are currently impractical to build because of 
their highly multiscale nature and because we assume the bars are perfectly rigid, and the 
junctions pin-like. Nevertheless, with the current rate of progress in rapid prototyping of 



complex structures, as exemplified by Kadic et al. (2012), it is hard to predict how well 



they might be approximated in the future. Also they show what is theoretically possible, 
and this should motivate the design of more realistic materials exhibiting a wide range 
of interesting behaviors. 

In the course of the analysis we construct two and three dimensional dilational ma- 



terials with an arbitrarily large flexibility window as defined by Sartbaeva et al. (2006) 
Dilational materials are perfect auxetic materials for which a dilation is the only easy 
mode of deformation: see the two-dimensional examples of figures Mjo) and|2]^c), and, 



in three dimensions, those of Biickmann et al. (2012) In my dilational material with an 



arbitrarily large flexibility window the cell size could (theoretically) expand from the size 
of a pea to the size of a house! 

Historically the study of what mechanisms can be achieved by linkages of rigid bars 
and pivots attracted considerable attention. In particular it was shown QKempe 1875 



lArtobolevskii 1964j |Kapovich and Millson 2002[) that any planar algebraic curve could 



be drawn by a linkage of bars and pivots (allowing for the bars to slide over one another 
so that the structure is not strictly planar). By contrast our focus is on the (macroscopic) 
mechanisms of periodic structures of bars and pivots. 
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Figure 2: Four examples of unimode materials. The example (a) of Larsen et al. (1997) 



which is a simplified version of a structure of Milton (1992), has an incremental Poisson 



ratio which is negative in the configuration shown. The examples (b) of Milton (1992) 



and (c) of Grima and Evans (2000) are dilators, for which the only mode of macroscopic 
deformation is a dilation. Example (d) is an expander which will be analyzed in detail 
in the next section. In (b), (c) and (d) the black regions are rigid polygons, which may 
be replaced by trusses of bars as shown in figure El The vertices in (a) where lines (rigid 
bars) meet, or in (b) where a line (rigid bar) meets a polygon, or in (c) and (d) where 
two or more polygons meet, are hinge joints (pivots). In each example an underlying 
Bravais lattice is marked by the vertices with red dots. 
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(a) 



Figure 3: As is well known a solid rigid polygon in two dimensions can be replaced by 
a truss of bars, as shown in (a) for a quadrilateral. Similarly a solid rigid polyhedron 
in three dimensions can be replaced by a truss of bars, as shown in (b) for a triangular 
prism: the basic idea is to subdivide the polyhedron into tetrahedra, and then replace 
the tetrahedra by bars along their edges. To simplify subsequent figures we will use rigid 
polygons (colored black) and rigid polyhedra (colored black and pink) rather than the 
equivalent truss. 

2 Realization of planar rectangular materials with 
arbitrary response 

In this section we only consider planar materials for which u and v are orthogonal and 
remain so as the material deforms. We call these rectangular materials even though 
the interior microstructure might not necessarily have rectangular symmetry. We let 
Ai = |u| and A2 = |v|. So for such materials the question becomes: what trajectories 
(Ai(£q), A 2 (to)) can be realized, as t varies in the interval t < t < t^. Such trajectories 
must lie in the quadrant Xi > and A2 > 0. We will show that every such trajectory 
can be approximately realized. Since we want to attach together unit cells which have 
different microstructures it is simplest to require that the parallelogram unit cell (which 
when periodically repeated forms the material) has all its microstructure remain inside 
the unit cell as the cell deforms and only touches the boundary of the unit cell at the 
vertices of the parallelogram which we call the support points. This is true for the 
microstructure of figure M^d) but not for the microstructures of figures Eta), E^b), and 
Et^c). This is not a restrictive assumption as for most microstructures, such as those in 
figures Ela), E^b) and Etc), we can take a large (finite) array of cells and attach bipod 
supports to the corners of the array and rescale the microstructure to obtain a new cell 
with almost the same properties as the original cell but for which the microstructure does 
not extend outside the unit cell: see figure E] to get the idea. 
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2.1 Some elementary operations on realizable functions 

To begin with we mainly concentrate on trajectories (Ai(to)> A 2 (io)) for which A2 is a 
single valued continuous function A2 = /(Ai) of Ai. As we will see there are many 
operations that can be done with realizable functions to obtain new realizable functions. 
We will not cover all of them, only those necessary to establish the main result of the 
paper, that all trajectories are approximately realizable. We start with the elementary 
(and obvious) ones, 

• Rescaling: If A 2 = /(Ai) is realizable then so is A2 = a~ 1 f(a\i), for any real 
constant a > 0, by rescaling the microstructure in the cell. 

• Rotation: If A 2 = f{X\) is realizable then so is Ai = /(A2) by rotating the cell 
(and its microstructure) by 90°. This is most useful when A 2 is a vertical distance 
which determines the horizontal distance Ai though the function /. 

• Reshaping: If A 2 = /(Ai) is realizable and / is continuous, A2 = rf(Xi) is re- 
alizable for any rational r = p/q, where p and q are integers. To see this, for 
simplicity in the case where p = 3 and q = 2, consider an array of np cells in the 
vertical direction and nq cells in the horizontal direction, as illustrated in figure H] 
with 71 — 4, where the bipod supports have leg length d (and have a pivot joint 
at the midpoint). By rescaling the microstructure of the material we can assume 
g = {nq)~ l f \nqh) . Also let the bipod supports scale as d — d /y/n, so they are 
small compared with Ai and A2 but large compared with the cell size. Then as 
n — > 00 we have nqh — > Ai and npg = rf{nqh) — > A 2 . So in the limit we see that 
A 2 = r/(A 1 ). 

• Addition: If A2 = /i(Ai) and A2 = /i(Ai) are both realizable and continuous, then 
A2 = /i(Ai)+/2(Ai) is realizable. Referring to figureEJ illustrated for the case n — 8, 
by rescaling the microstructure of the material we can assume g\ = (n)~ l fi{nh) 
and gi = (n)~ 1 f2(nh). We again let the supports scale as d — do/y/n. Then as 
n —7- 00 we have nh — > Ai and ng\ + ng 2 = fi(nh) + fzinh) — > ^2- I n the limit we 
obtain A 2 = /i(Ai) + / 2 (Ai). 

2.2 Some elementary realizable functions 

We obviously need some elementary realizable functions to begin building more compli- 
cated ones. Let us examine in more detail the response of the expander of figure EJ^c). 
The unit cell is shown in more detail in figure [6(a), which defines various parameters. 
Elementary trigonometry shows that 

A x = 2a sin a + 2e cos a, A 2 = 2acosa, (2.1) 

implying that (Ai, A2) lies on the portion of the ellipse 

(aAi - e\ 2 ) 2 = a 2 (Aa 2 - \\), (2.2) 

with 2y/e 2 + a? > Ai > 2s and 2a > A2 > 2ae/\/e 2 + a 2 where the last constraint ensures 
that the triangles remain in the cell. 
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Figure 4: The configuation for reshaping a cell, showing that for any rational r = p/q, 
A2 = (p/o)f(^i) is realizable if A2 = f(Xi) is realizable. The configuration is illustrated 
for the case n = 4, p = 3 and q = 2, and one should take the limit n — )■ 00. The leg 
length d of the bipod corner supports has to be scaled appropriately as discussed in the 
text. 
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Figure 5: The configuration for an adder, combining two materials (red and blue) to 
obtain a metamaterial with a response function which is the sum of the response functions 
of the two materials. The bipod corner supports have leg length d which is large compared 
with h, gi and #2 yet small compared with Ai and A2. 
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Figure 6: Shown in (a) is the basic geometry of a cell of an expander. Here e should 
be chosen very small. The black right angled triangles are rigid. Shown in (b) is an 
expander of n — 10 cells, with the legs of the corner triangles chosen slightly longer, so 
they can be used as supports. 
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In particular let us take n of these elements, arranged as in figure HD^b), to form a 
new cell with dimensions (A'^Af,) and let us set e = eo/n 2 and d = d /n 2 keeping a 
independent of n. Then given a domain D of operation, Aj~ < X[ < Xf, we choose 
n sufficiently large that the minimal extension of the expander is less than Aj~ and its 
maximum extension is greater than A J. Then for any given \[ in D we have 

< (aAi — e\ 2 ) < aAi < a\[/n < a\f /n, A' 2 < A 2 < A' 2 + d /n 2 . (2.3) 

These estimates in conjunction with (12 .2\\ imply that for any X[ in D, A' 2 converges 
to 2a as n — > oo, independent of X[. In other words any positive constant function 
is (approximately) realizable on a given domain D. Such a cell behaves like an ideal 
expander: when n is very large we can adjust the horizontal length as we please (within 
the domain D) and the vertical length remains (essentially) unchanged. From a practical 
viewpoint we caution that the structure will likely be very sensitive to small imperfections 
when n is large. 

Another important cell is a dilator. An example of one, based on the chiral geometries 



of |Milton (1992)| |Prall and Lakes (1996)1 and |Mitschke et al. (2011)[ is shown in figure 
[7J in four states of deformation. If the innermost square has side length £ and the right 
angled triangle has height h and base h < £, then the cell dimensions are 



\ 1 = A 2 = t = \/2h 2 + £ 2 - 2V2h£cos(6 + tt/4), (2.4) 
having a minimum value of 

to = V2/i 2 + t 2 - 2h£ (2.5) 
at 9 = and a maximum value of 

t+ = £ + y/2h (2.6) 

at 9 — 37r/4 corresponding to case (b) in figure [71 We call (t^t^) the performance 
range of the dilator. Rescaling the microstructure allows us to get performance ranges 
(fctg, Hq) for any constant k > 0. 

Now by combining four expanders with this basic dilator cell, as sketched in figure El 
one can obtain dilators having an arbitrarily large performance ratio / tg of maximum 
extension to minimum extension. To see this, suppose each expander has n unit cells and 
let the dimensions of the inner basic dilator scale as £ = £o/n, h = h /n and b = bo/n, 
let the dimensions of the expander scale as a = a /n, e = eo/n 2 , and let the dimensions 
of the bipod supports scale as d = do/n. Choose a in the range 



2h 2 + £ 2 - 2h £ /V2 <a < (£ + V2h )/2, (2.7) 

so that the inner basic dilator can never reach its maximum extension, but the expander 
can attain its minimal extension at a = and can have an angle a at least as large as 
45°. (This inequality implies £o and ho must be chosen with 6ho£o > 2h\ + £\ as may 
be satisfied with ho = £o, for example). Choose d > \/2ao so that the microstructure 
is always contained in the cell. When the expanders are fully compressed they have a 
length of 2ne = 2eo/n, and since the maximum possible dimensions of the inner basic 
dilator scale as 1/n and the supports scale as 1/n it follows that the minimum extension 
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scales as 1/n. By contrast with a = 45° the expanders have an extension of \/2ao in 
the limit n — > oo. Therefore the performance ratio of maximum extension to minimum 
extension can be as large as one pleases with this microstructure by choosing n large 
enough. Thus the function A2 = Ai is realizable, and by addition with an ideal expander 
so is the function A 2 = Ai + c for any constant c > 0. 




Figure 7: Four successive configurations of a basic dilator cell. The maximum area of 
the dilator enclosed by the square formed by the outermost points occurs in (b) where 
this square is aligned with the innermost square 

Now that we have dilators we can build many things. 

• To build a material where in a domain A2 depends linearly on Ai with negative 
slope consider the microstructure of figure |9j In the cell there are n dilator cells in 
the horizontal direction and n dilator cells in the vertical direction. We only show 
one corner of the cell as the microstructure in the other corners are just reflections 
of this microstructure about the x\ and £2 axes. We scale the geometry with 
a = ao/y/n, b = a /y/n, c = c /\/n, d = d /\/n, and (ft = (f)o/\/n while varying 
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9 so that 6 + is independent of n. Since by trigonometry a + 7 + = 7r and 
# + + 7 — /3 = 7rwe deduce that a + f3 = <fi = (fio/y/n. So as one might expect, the 
angles a and (3 must go to zero as n — )■ 00. We can let a = ao/y/n and /3 = Po/V™ 
where ao > and /3o > depend on the deformation and ao + 0o = 4>o- Then the 
length nh of the dilator segment in the horizontal direction satisfies 

nh = 2nb sin a = 2 y/nbo sin (cto/ — >■ 26o«o as n -> 00, (2.8) 

while the length ng of the dilator segment in the vertical direction satisfies 

ng = 2ncsin/3 = 2 y/nc sin (/3 /y/n) — > 2c /3o as n — > 00. (2.9) 

We restrict the range of operation so Ai > rj and A2 > 1] for some small rj > 0. 
Then (Ai, A2) must approach (nh,ng) since the corner cell geometry scales roughly 
as l/i/n. Finally since ao + Po = <po we conclude that, in the limit n — > 00, 

A 2 = 2c o 0o - (cd/6 )A 1 . (2.10) 

Since the ratio Co/60 can be any positive value we please, it is clear that we can get 
any negative slope. By rescaling it is clear that the constant can take any desired 
positive value. 

To build a material where in a domain D, < A^f < Ai < A^, A 2 depends quadrat- 
ically on Ai consider the microstructure of figure [10J In the cell there are n 2 dilator 
cells in the horizontal direction and n 3 dilator cells in the vertical direction. Like 
in the previous example we only show one corner of the cell as the microstructure 
in the other corners are just reflections of this microstructure about the x\ and X2 
axes. Since a and b are sides of a triangle with opposite angles and 7r/2 — a the 
law of sines implies sin = (a/b) cos a. To ensure that 7 = when a = 0, in which 
case (3 = 6, we choose sin# = a/b. Then we have 

g = 2c sin 7 = 2csin(6 l — 0) 

= 2c 



sin^y 1 — (a/b) 2 cos 2 a — (a/b) cos 6 cos a 

(2ca/b) Ul - (a/b) 2 + (h/2b) 2 - - (a/b) 2 )(l - (h/2af 
eh 2 

+ 0(h 4 /a A ). (2.11) 



Aay/b 2 - a 2 
Now we take a scaling with 

a = ao/n, b = bo/n, c = co/n, e = Eo/n, d = do/n, (2-12) 

while keeping 9 independent of n. Since the corner microstructure shrinks to zero 
as n — > 00, given Ai in the domain D, n 2 h must converge to Ai. Hence h/a = fin/a® 
is of the order of 1/n and converges to zero as n — > 00. So from (12.1 ip we see that 
n 3 g approaches c (Ai) 2 /[4a A/&o — a o] (with the leading correction terms in (12. lip 
contributing a term of order 1/n), which in this limit can be identified with A 2 
because the corner microstructure shrinks to zero. Hence we deduce that 



\ 2 = k\l, with k = Co/[4a Jbl-a 2 l }. (2.13) 
By rescaling it is clear that the constant k can be chosen to have any desired value. 
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b 

Figure 9: The corner of a cell which realizes a linear dependence A2 = r — sAi (where r > 
and s > 0) to an arbitrarily high degree of approximation within a domain Ai > f] > 
and A2 > 77 > . Here the green cells are dilator cells. The insert shows the geometry of 
the rigid quadrilateral, defining the angles <fi and 9 
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Figure 10: Design of the squarer. The corner of a cell which, within a domain, realizes 
a quadratic dependence A 2 = k\\ to an arbitrarily high degree of approximation within 
a domain. Here the green cells are dilator cells. The insert shows the geometry of the 
rigid triangles, defining the included angle 9. 
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• A variant of the squarer microstructure as illustrated in figure [TTJ is obtained by 
setting a = b and 6 = n/2. We take n 2 dilator cells in the horizontal direction and 
n 2 dilator cells in the vertical direction. In this geometry we have h = 2a sin a and 
g = 2c sin a. This implies that for Ai in a domain D, < A^f < Ai < Xf, in the 
n — > oo limit we have 

A 2 = fcAi, with k = c/a. (2.14) 
Of course for rational k we could have obtained the same response by reshaping. 




Figure 11: The corner of a cell which, within a domain, realizes a linear dependence 
A2 = k\i to an arbitrarily high degree of approximation within a domain. Here the green 
cells are dilator cells. The insert shows the geometry of the right angled rigid triangles. 

2.3 Some additional operations on realizable functions 

Let us define T> as the set of all finite connected subintervals which are contained strictly 
inside the positive real axis, which may be open or closed at either end: 

V = {(a, 6), (a, b], [a, 6), [a, b] \ 0<a<b< 00}. (2.15) 

We have the following two additional operations on realizable functions, 

• Subtraction: Not only can we realize the sum of two realizable continuous 
functions /i(Ai) and /2(Ai), defined (and operational) on an open domain D G V, 
but we can also realize the difference f 3(^X1) = fi(Xi) — f 2(^X1) on any subinterval 
of D provided /3(D) G V. We emphasize that /i(Ai) — /2(Ai) could be negative 
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outside D. To see this realization consider the microstructure of figure [121 In the 
cell there are n dilator cells in the horizontal direction and n dilator cells in the 
vertical direction. We only show one corner of the cell as the microstructure in 
the other corners are just reflections of this microstructure about the x\ and x 2 
axes. Let I = [Ar, A^] be a closed subinterval of D . We choose to scale the bipod 
support leg length d as d = d Q /\/n. The inequalities 

nh < Ai < nh + 2d = nh + 2d /Vn (2.16) 

imply as n — > oo that nh converges to Ai from below. Also it implies there exists 
a no (which can made explicit) such that nh G D for all Ai G / when n > n$. 
Then f\{nh) and f2(nh) are defined. By rescaling and by a 90° rotation of the 
microstructure of the material we can assume s = {n)' 1 f\{nh) and r = {n)~ l f 2 (nh) . 
By construction g = s — r = (n)^ 1 f 3 (nh) . So from the inequalities 

ng < A 2 < ng + 2d = ng + 2d /\fn ( 2 -17) 

we see that ng = fz{nh) must converge to A 2 in the limit n — > oo. Therefore, since 
/3 is continuous on D, in this limit A 2 = /s(Ai) for Ai G /. 

• Composition: We can also compose two realizable continuous functions / 2 (A 2 ) 
and /i(Ai) to obtain a realizable function A3 = /s(Ai) = / 2 (/i(Ai). We are given 
open domains Di G T> and D 2 G T>, such that f\ and / 2 are defined on D\ and 
D 2 , respectively, and /i(-Di) C D 2 . Also we assume fs(Di) G V. We consider a 
subinterval I = [A^, A^] C T)\. The composition is done through the microstructure 
of figure [T3j In the horizontal direction there are n dilator cells and in the vertical 
direction there are n cells of rescaled material corresponding to / 2 . At the junction 
there is the rescaled material corresponding to fx rotated by 90°. We choose to 
scale the bipod support leg length d as d — d /y/n, so the inequalities (12. 16ft hold 
and imply there exists a n$ such that nh G D\ for all Ai G I when n > no- 
Then fi(nh) and fz{nh) are defined. By a 90° rotation and rescaling we can realize 
t = (n)^ 1 fi(nh) and by rescaling realize g = (n)~ 1 / 2 (nt) = fs(nh). The inequalities 

ng < A 2 < ng + 2d = ng + 2d /y/n (2.18) 

imply that ng = fz{nh) must converge to A3 in the limit n — > 00. Therefore, since 
fs is continuous on D, in this limit A 2 = feiXi) for Ai G I. 



2.4 Realizability of arbitrary functions 

We first prove that one can on any subinterval I = [A^f, A^] of a fixed domain D G T> 
0<A ] ^<Ai<A ] f realize any polynomial 

A 2 = p(Ai) = a + aiAi + a 2 A? + a 3 Ai + . . . + a n A™, (2.19) 

which is bounded below by a positive constant on this domain. We prove this assertion 
by induction. Suppose it is true for polynomials p of degree n less than or equal to 
2m, where m > 0. This is certainly true when m = 0. Now to avoid carrying around 
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Figure 12: The configuration for the corner of a subtractor, combining two materials 
(blue and red) to obtain a metamaterial with a response function which is the difference 
of the response functions of the two materials, within the domain of operation. Here the 
green cells are dilator cells. 
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Figure 13: The configuration for a composer, combining two materials (blue and red) to 
obtain a metamaterial with a response function which is the composition of the response 
functions of the two materials. Here the green cells are dilator cells. 
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superfluous constants let us suppose units of length (say centimeters) have been chosen, 
so for example A2 = 1 means A2 equals one centimeter. Similarly A2 = A 2 means this 
relation holds when both Ai and A2 are measured in centimeters, i.e. A2 = kX\ where 
k = 1cm -1 . Since we can realize A™ +1 , Ai + 1, and (by the construction of the squarer 
of figure [TOl) A 2 it follows by composition that we can realize (even when m = 0) 

A 2m + 2 = ( A m+1)2 ^ + ^2 m+ 2 = y^n+2 + ^ + 2 )\^ + gf(Al), (2.20) 

where by the binomial theorem g(\i) is a polynomial of degree 2m. By rescaling we can 
also realize &iA 2m+2 and &2(Ai + l) 2m+2 for any positive constants b\ and 62- 

Now suppose we are given any polynomial q{X\) of degree 2m + 2 or less. Using (I2.20p 
we can express it in the form 

g(A0 = Cl \l m+2 + c 2 (Ax + l) 2m + 2 + r(Ai), (2.21) 

where r(Ai) is of degree 2m or less, and the constants C\ and Ci could have either sign, 
or be zero. Now there exists a sufficiently large constant c > such that c + CiA 2m+2 , 
c + C2(Ai + i) 2m + 2 and c + r(Ai) are bounded below by a positive constant on the interval 
/. Hence using the fact that sums or differences of realizable functions are realizable 
(provided the difference is bounded below by a positive constant on I) we deduce that 
each of these three functions is realizable on J, and hence too is their sum s(Ai) in terms 
of which 

g(Ai) = s(Ai) -3c. (2.22) 

Again this is the difference of two realizable functions, and will be realizable if q(Xi) is 
bounded below by a positive constant on the domain /. This proves the assertion is true 
for all polynomials of degree n < 2m + 2, and by induction for polynomials of any degree. 



Now by the Weierstrass approximation theorem (see, for example, Keener (2000) ) 
given any continuous function f(\\) on J, we can approximate /(Ai) arbitrarily closely 
by polynomials. Specifically, given any e > there is a polynomial p(\\) so that 

max|/(A 1 )-p(Ai)|<e. (2.23) 
\iei 

Since we can realize p(\i) arbitrarily closely it follows that we can realize /(Ai) arbitrarily 
closely. 

2.5 Realizability of any trajectory 

Let R(Xi, Aj, X2) denote the rectangle 

< Ar < Ai < Xi, < A 2 < A 2 < Xt, (2.24) 

and let 1Z denote the set of all such rectangles. Here we show that any continuous 
trajectory (Ai,A2) = (/i(^o), ^2(^0)) taking values in a rectangle in 1Z is realizable. We 
assume the parameterization has been chosen so to has dimensions of length, and to 
increases from t$ > to £q > £q along the trajectory, and that fi(t ) and ^2(^0) are 
defined, bounded, continuous and positive on the closed interval / = [io>*o1- ^ n ^ s 
implies fi(I) € T> and /2OO £ D. The geometry which accomplishes the realizability is 
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that in figure UM The blue material, of which there are n cells in the horizontal direction, 
is a rotation by 90° of the material which realizes /i(io) rescaled so h = n^f^nt). The 
red material, of which there are n cells in the vertical direction, is the material which 
realizes / 2 (to) rescaled so g = n~ 1 f 2 (nt). We choose t = t /n and we scale the bipod 
support leg length d as d — d / \fn. Then for t G / we have 

nh < Ai < nh + 2d = nh + 2d Q /yfn, (2.25) 

which implies f\(to) = nh converges to Ai as n — > 00. Similarly /2(^o) = n 9 converges 
to A 2 as n — )■ 00. Therefore the trajectory (Ai,A 2 ) = (/i(^o) ; /2(^o)) is realizable. An 
interesting corollary is that the trajectory can self-intersect. At points of intersection 
knowledge of Ai and A 2 is not sufficient to determine t Q and hence to determine what 
subsequent deformations are possible: one needs to keep track of the hidden variable to- 




Figure 14: Given two planar rectangular materials (blue and red) having arbitrary re- 
sponse functions one can construct an rectangular metamaterial having any desired re- 
sponse trajectory (Ai(t), A 2 (t)) within a desired domain. Here the green cells are dilator 
cells. 



3 Realizability of oblique materials with an arbitrary 
response 

Let us introduce polar coordinates (r,9) such that (Ai,A 2 ) = (r cos 9, r sin 6). Let 
A{r~ , r + , v) denote the sector 

<r -<r<r + , < v < 9 < n/2 - u, (3.1) 
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and let A denote the set of all such sectors. Since, with rectangular materials, we are 
free to realize any trajectory (Ai(i), A 2 (£)) taking values in a rectangle in TZ, we can in 
particular realize any trajectory where r increases strictly monotonically with t such as 
the blue trajectory in figure [J5j In other words we can realize any continuous desired 
trajectory 9 = 0(r) taking values in a sector in A, with Ai = e\{r) and A2 = e 2 (r) for 
some continuous functions &i{r) and e 2 (r). 




Figure 15: A desired function 9(r) taking values in a domain < 9~ < 9 < 9 + < n/2, 
< t~ < t < t + corresponds to a trajectory (Ai(t), A 2 (t)) (for example the blue one 
shown here) where (r, 9) are the polar coordinates of (Ai, A 2 ). 

Next consider the geometry of figure [16] which we call an angle adjuster. The yellow 
cells are rectangular material. Each of the four arms in the cross in is 2 cells wide and 
n cells long (in the figure n = 7) and the rectangle at the center of the cross is two cells 
by two cells. Here by rescaling h = (2n) -1 ei(2w), g = (2n)~ 1 e 2 (2nr) and 9 = 0(2nr) 
in which r is the length of the cell diagonal and tan# = g/h, as shown in the insert. 
Thus the two horizontal arms have lengths e\(2nr)/2 and the vertical arms have lengths 
e 2 (2nr)/2. We choose to scale the bipod support leg length d as d = d§j \Jn. Thus in the 
limit as n — >■ 00 the width of the arms shrinks to zero, and the size of the bipod supports 
shrink to zero, and 2nr — > t while a — > 29. So, by the continuity of 0, in the limit we 
have that a = 2<p(t). Note that by rescaling of the microstructure we can also realize 
a = 2<p(ct) for any constant c > 0. 

With this angle adjuster we can now obtain oblique materials (i.e. periodic materials 
with an underlying oblique lattice) with an arbitrary response through the microstructure 
of figure [171 Now u and v are no longer orthogonal, and we let r, with < r < 7r denote 
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Figure 16: By taking a material (shown yellow here) with an rectangular response one 
can obtain an angle adjuster with a cell which is rhombus. The insert shows the geometry 
of one cell. 
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the angle between these vectors, which varies as the cell deforms. In the corners is one 
cell (yellow) which is an angle adjuster. On the top and bottom there are n cells of blue 
rectangular material and on the sides there are n cells of red rectangular material. We 
define B as the set of all finite connected subintervals contained strictly in [0,7r], which 
may be open or closed at either end: 

B = {(a, b), (a, b], [a, b),[a,b] \0<a<b< vr}. (3.2) 

Suppose we are given any two functions /i(t ) and /^o) which are continuous on an 
open interval D 6 V, with fi(D) 6 V and /2OD) £ T>, and a function (fi(t ) which is 
continuous on D with <p(D) G B. We can, by rescaling realize a blue rectangular cell with 
h = n~ l fiint), a red rectangular cell with h = n~ 1 f 2 (nt) and a yellow angle adjuster cell 
with a = <p{nt) for nt in any closed subinterval / = [t$ , t$ } C D. We choose t = to/n and 
scale the bipod support leg length d as d = d /y/n. Then in the limit n — > 00 we have 
Ai = |u| — > nh, A2 = |v| ng and a — > r. In this limit it then follows that Ai = f\{to), 

= ^2(^0) an d t = 4>(t ) for all t I- We have thus realized any desired response. 

Finally, to make sure the trajectory does not extend outside [toi^ol we choose the 
dilators at the corners in figure [T71 so that their performance range is exactly nt e (t^, t$ ) 
or, alternatively, given the dilator performance range we parameterize the trajectory to 
match it. The endpoints of the trajectory t = t$ and to = ^o" m ight then not be 
realizable. 




Figure 17: By combining two rectangular materials (blue and red) and a rhombus shaped 
angle adjuster (yellow) one can obtain a oblique material having an arbitrary response. 
The green cells near the corners are dilator cells. 
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4 Realizability of three-dimensional orthorhombic ma- 
terials with an arbitrary response 



A key element in the construction of three dimensional materials with any desired re- 
sponse is the construction of panels from two dimensional rectangular or parallelogram 
cells lying in the x\ — Xi plane. An example of the construction of a panel from the 
expander cell is given in figure ITST a). The first step is to extend the microstructure a 
short distance r (short compared to the cell size) in the positive £3 direction, so rigid 
lines become rigid faces, rigid triangles become rigid triangular prisms, rectangles become 
rectangular prisms, and so forth. Pivot junctions become edge junctions. The next step 
is to modify the structure so it tapers to a point at the four supports. The final step is 
the bend the tips of these supports a short distance s in the negative x 3 direction, so that 
the microstructure is strictly contained in a box having the four supports as corners of 
one face. The final construction may have rigid flat rectangles as in [151(b). If one a desires 
a construction that only involves rigid bars then one should first make the rectangle into 
a thin triangular prism as shown in figure [TBT c) and then make the transformation to a 
truss of bars as in figure Mp). 



(a) 



X, 





t 



(c) 




Figure 18: As shown in (a) two dimensional cells can be extended into the x 3 -direction 
to obtain panels, which can be used to construct three dimensional cells. The panels are 
tapered at their four supports, and the tapered ends are bent slightly in the x 3 -direction, 
so that all the microstructure is to one side of the plane through the four support points. 
To avoid confusion those faces of the rigid objects parallel to the X3-axis are colored in 
pink. Any rigid flat rectangles in the final construction as in (b) can be transformed into 
thin triangular prisms as in (c) and then replaced by a truss of bars if desired. 



Following the construction of Buckmann et al. (2012) and as shown in figure IT91 one 



can assemble 5 (or 6) panels based on the cell geometry of figure [7] to form a cubic cell 
which acts as a three-dimensional dilator. To obtain three-dimensional dilator cells which 
allow arbitrary expansion one could use the cell structure of figure [5] as the basis for the 
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panels. A periodic dilational material of these cubic cells would have an arbitrarily large 
flexibility window as defined by Sartbaeva et al. (2006) 



(a) 




(b) 




Figure 19: Shown in (a) is the panel which corresponds to the two-dimensional dilator 
cell of figure [7J To avoid confusion those faces of the rigid objects parallel to the a^-axis 



are colored in pink. Following Buckmann et al. (2012) copies of these can be assembled 



as the faces of a cube to form a three-dimensional dilator cell, as illustrated in (b). Only 
five, rather than six, panels are shown. In fact these are all that are needed. For clarity 
two of these panels are colored in brown and the three dimensional structure of each 
panel is omitted. They are oriented so the microstructure lies strictly inside the cube 
formed by the eight support points. 

Now we can use these dilators to construct three-dimensional orthorhombic materials 
with an arbitrary response. By orthorhombic we mean that the vectors u, v and w are 
mutually orthogonal and remain so as the material deforms. Let the cell of the material, 
that we will construct, have sides of lengths Ai = |u|, A2 = |v| and A3 = |w| in the Xi, 
x 2 , and £3 directions. Let P(X[, Xf, A^, A^, A3 , A3 ) denote the rectangular prism 



< A^ < Ai < A^ 



< A 2 < A 2 < \J, 



< A, < A3 < A 



3 • 



(4.1) 



and let V denote the set of all such rectangular prisms. Here we show that any trajectory 
(Ai, A 2 , A 3 ) = (fi(t Q ), f2(t ), ^3(^0)) taking values in a rectangular prism in V is realizable. 
We assume the parameterization has been chosen so to has dimensions of length, and to 
increases from t$ > to t$ > t$ along the trajectory, and that /i(to), ^2(^0), and 
^3(^0) are defined, bounded, continuous and positive on the closed interval / = [t^ ,£q]. 
This implies f\{I) G T>, /^(-O G T> and /3(/) G T>. The corner of the geometry which 
accomplishes the realizability is shown in figure |2DJ This geometry is a generalization of 
the two dimensional geometry of figure dH In the x±, x 2 , and £3 directions there are blue, 
red, and cyan tubes, respectively, which are each n units long, where each unit consists of 
four panels joined to make a square tube. Each tube retains its square cross section due 
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to its contact with the green dilator cells. Choose t = to/n, where to G I, and scale the 
tripod support leg length d as d = d^j \fn. From the results of section 2.4 and rescaling 
we can realize h = n~ 1 fi(nt), g = 71/2 and s = n" 1 f$(nt). Then for to G / we have 

nh < Ai < nh + 2d = nh + 2d /y/n, (4.2) 

which implies /i(to) = nh A i as n — >■ oo. Similarly it follows that ^2(^0) = ng X2 and 
/3(*o) = ™s -> A 3 as n ->> 00. Therefore the trajectory (Ai, A 2 , A 3 ) = (fi{t ), f 2 (t ), /3^o)) 
is realizable. 




Figure 20: Corner of the cell which realizes any desired orthorhombic response. The 
green cell is a three-dimensional dilator cell and the blue, red, and cyan tubes, with 
square cross section, have panels based on cells with the needed response. 



5 Realizability of three-dimensional triclinic materi- 
als with an arbitrary response 

The construction of three-dimensional triclinic materials with an arbitrary response, is 
similar in some respects to the two-dimensional construction of figure [T71 We let <p\ 
denote the angle between u and v, 02 the angle between v and w, and 03 the angle 
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between w and u. Assuming the vectors u, v and w are not coplanar (since otherwise 
A = F T F is singular) these angles satisfy 

< 0i < 7T , < 2 < 7T, < 03 < 7T, 0i + 02 + 03 < 27T. (5.1) 

The first step is the construction of a three-dimensional angle adjuster as shown in figure 
l2~TT a). just before the points Di and D2 are joined. The construction of one corner of the 
cell then proceeds as shown in figure l2TT b) . See the caption for more details. In the xi, 
X2, and 23 directions there are blue, red, and cyan tubes, respectively, which are each n 
units long, where each unit consists of four panels joined to make a square tube. Three 
additional rhombi, with side length mt (composed of m x m orange, m x m magenta, 
and m x m gold panels) can be added to the corner of the cell of figure 121T b) to form 
a rhombohedron, which can be copied 7 times and translated to the other 7 corners of 
the cell and joined to the tubes with dilator cells, in similar way to the construction of 
l2~iT b). Now if one has two line segments of length nh + 2(m — l)t (the length of the 
interval between C and the equivalent point in the vertical direction at the top corner of 
the rhombohedron above it) joined to two line segments of length ng + 2(m — l)t where 
the angles between line segments are either a or 7r — a, then those line segments must 
form a parallelogram (assuming, as is the case here, that "opposite angles are the same" 
to avoid the case of a warped isosceles trapezoid). The corresponding line segments from 
all sides form a parallelepiped, with angles of a, tc — a, (3, tc — (3, 7, and tt — 7, with 
side lengths of nh + 2(m — l)t, ng + 2(m — l)t, and ns + 2(m — l)t. A distance t from 
the corners are the attachments to the tripods, having leg length d. So the support is at 
most a distance d from these attachments. Thus we have the inequalities 

nh + 2(m - 2)t - d < Ai < nh + 2(m - 2)t + d, 
ng + 2(m - 2)t - d < A 2 < ng + 2(m - 2)t + d, 
ns + 2(m - 2)t - d < A 3 < ns + 2(m - 2)t + d, (5.2) 

where Ai = |u|, A 2 = |v|, and A 3 = |w| are the lengths of the sides of that parallelepiped 
which has the eight support points (each at the end of a tripod) as vertices: see figure 

m 

Next suppose we are given any three functions /i(£q), ^2(^0) an d ^3(^0) which are 
continuous on an open interval D 6D, with fi(D) G T> for i = 1,2, 3, and three functions 
0i (to), 02 (£0) an d 03 (to) which are continuous on D E V, with (f)i(D) G B for % = 1,2,3, 
and additionally with ip(D) G B where 

<p(t Q ) = (0i (t ) + 2 (t o ) + 3 (£o))/2. (5.3) 
We can realize, by rescaling, 

• a blue rectangular panel with h = ra _1 /i(nt), 

• a red rectangular panel with g = n~ 1 f2{nt) 

• a cyan rectangular panel with s = n~ 1 f2(nt) 

• a orange angle adjuster panel with a = 0i(nt) 
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• a magenta angle adjuster panel with (3 = faint) 

• a gold angle adjuster panel with 7 = faint) 

for nt in any closed subinterval I = [to,to] C D. We choose t = to/n and scale the 
tripod support leg length d as d — do/y/n. Since faiD) E B for i = 1,2,3 it follows that 
there exists an angle n such that 

n < fait ) < n - 77, i = l,2,3. (5.4) 

Thus, for to £ I none of the angles a, ix — a, /3, ir — (3, 7, and tt — 7 are less than 17. We 
choose m independent of n and such that 

m > 1 + l/tan(r//2), (5.5) 

to ensure that none of the green dilation cells touch one another (and possibly jam the 
structure) as the cell deforms: see figure [23j Then in the limit n — > 00 (15.2ft implies 
Ai nh, A2 — > ng and A3 — > ns. Likewise, because the size of the supports shrinks to 
zero, 7~i — > a, r 2 — > f3 and r 3 — > 7. This, implies that in the limit n — > 00, 

Ai = /i(io), Ti = fait ), i = l,2,3. (5.6) 

We have thus realized any desired triclinic response. Finally, to make sure the trajec- 
tory does not extend outside [£q ,£q] we choose the dilators in figure [2lTb) so that their 
(rescaled) performance range is exactly nt 6 (£q ,t|j~) or, alternatively, given the dilator 
performance range we parameterize the trajectory to match it. The endpoints of the 
trajectory to = £q and to = to might then not be realizable. If an actuator is placed 
within one of these dilator cells in the material then it should be possible to adjust to to 
any desired position along the trajectory. 
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Figure 22: The vertices of this parallelepiped are the support points, each at the end of 
a tripod. The microstructure inside the cell, including the tripods, is not shown, but can 
be visualized using figure |21j that corner corresponds to the corner here where the three 
dashed lines meet, although the angles are slightly different due to the presence of the 
tripods. The angles 0i, 02 and 03 are the actual angles between lines (not the angles 
between their projections) and thus sum to less than 360°. When they are all 90° the 
parallelepiped becomes a rectangular prism. 




Figure 23: The worst configuation of the green dilator cells relative to an orange angle 
adjuster. They avoid collision if tan(a/2) > \/{m — 1). Hence follows the condition (15. 5ft . 
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